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Abstract. Let C C be a subfield of the complex numbers, and let D 
be the ring of K-\\aeai differential operators on R = K[x^... ,x„]. If M 
and A*^ are holonomic left D-modules we present an algorithm that computes 
explicit generators for the finite dimensional vector space }lom£){M, N). This 
enables us to answer algorithmically whether two given holonomic modules 
are isomorphic. More generally, our algorithm can be used to get explicit 
generators for Ext^(M, Af) for any i. 



1. Introduction 

Let D = Dn = K{xi, . . . , di, . . . , 9„) denote the n-th Weyl algebra over a 
computable subfield K C C, i.e. elements of K can be represented with a finite set 
of data, their sums, products and quotients can be calculated in a finite number of 
steps, and there is a finite procedure that determines whether a given expression 
of elements of K is zero or not. Let Hom£)(M, TV) denote the set of left ZJ-module 
maps between two left D-modules M and N. Then HomD(M, A^) is a iiT-vector 
space and can also be regarded as the solutions of M inside TV in the following 
way: Given a presentation M ~ 0"^° /D • {Li, . . . , Lr^}, let S denote the system of 
vector- valued linear partial differential equations, 

5 = {Li . / = . • • = L,, . / = 0}, 

and let So^S*; iV) denote the A^-valued solutions / £ N''° to S. Then the ho- 
momorphism space Homu (1?''° /_D • {Li, . . . , Ln}, iV) is isomorphic to the solution 
space 801(5*; N) where the identification is as follows. A homomorphism ip in 
E.omD{D^° / D ■ {Li, . . . , Lr^}, N) corresponds to the solution [ip{ei), . . . , (p{ero)]'^ £ 
TV" of 5, while a solution / = [/i, . . . , frg]'^ G of S corresponds to the homo- 
morphism which sends et to fi. 

If M and iV are holonomic, then the set IIom£)(M, N) as well as the higher de- 
rived functors Ext]j{M, N) are finite-dimensional i^T- vector spaces. In this paper, 
we give algorithms that compute explicit bases for IIomD(M, N) and Ext^(M, N) 
in this situation. Our algorithms are a refinement of algorithms given in p2| , which 
were designed to compute the dimensions of }iomu{M, N) and Ext^(M, A'') over 
K. Algebraically, our problem of computing a basis of homomorphisms is easy to 
describe. Namely, since a map of left £>-modules from M to N is uniquely deter- 
mined by the images of a set of generators of M, we must simply determine which 
sets of elements of N constitute legal choices for the images of a homomorphism (of 
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a fixed set of generators of M). It is perhaps surprising tliat this is a difficult com- 
putation. One of the reasons is that Home {M, N) lacks any Z?- module structure 
in general and is just a _fC- vector space. 

In recent years, one of the fundamental advances in computational _D-modules 
has been the development of algorithms by Oaku and Takayama ^, |lOj to compute 
the derived restriction modules Torf {D/{xi, . . . , Xd} ■ D, M) and derived integra- 
tion modules Torf {D/{di, . . . , dd] ■ D, M) of a holonomic Z?-module M to a linear 
subspace xi = ■ ■ ■ — Xd — 0- We give a summary of these algorithms in the ap- 
pendix. These algorithms have been the basis for local cohomology and de Rham 
cohomology algorithms |l^ and have been extended to algorithms for derived 
restriction and integration of complexes with holonomic cohomology by the second 
author [|6|. 

Similarly, the algorithm of to compute the dimensions of IIom£)(Af, N) and 
Ext^(M, iV) is also based on restriction by using isomorphisms of Kashiwara and 
Bjork 1^. These isomorphisms are, 

(1) Extl,(Af, N) - Tor^_,(Ext2,(Af, D),N), 
which turns an Ext computation into a Tor computation and 

(2) Torf (M', N) ~ Torf " (i?2„/{x, -y^,^, + *,}^Li ■ ^2„, t(M') M N), 

which turns a Tor computation into a twisted restriction computation in twice 
as many variables (an explanation of the notation used above can be found in 
Section ||). 

In this paper, we will obtain an algorithm for computing an explicit basis of 
Ext^ (M, N) by analyzing the isomorphisms (|l|) and (|^) and making them com- 
patible with the restriction algorithm. In Section 2, we present a proof of isomor- 
phism (|l|) adapted from In Section 3, we give an algorithm for computing 
Honi£i(Af, iV) in the case N = K[xi, . . . ,a;„], which is used to compute polyno- 
mial solutions of a system S. In Section 4, we give an algorithm for the case 
N = K[xi, . . . , a;„][/~^], which can be used to compute rational solutions of S. In 
Section 5, we give our main result, which is an algorithm to compute IIom£)(Af, N) 
for general holonomic modules M, N. In Section 6, we give a companion algo- 
rithm which computes Ext^ (M, N) and their representation in terms of Yoneda 
Ext groups. In Section 7, we give an algorithm to determine whether M and N are 
isomorphic and if so to find an isomorphism. We also give an algorithm to com- 
pute the endomorphism ring End£)(A/), the algebraic group 1sod{M), and some of 
their basic properties. In the appendix, we review the restriction and integration 
algorithms. The reader may refer to Algorithm |8.6| for a d iscussion of the restric- 



tion algorithm and the t^-filtration, and to Algorithm |8.8| for a discussion of the 
integration algorithm and the ^-filtration. Finally, the algorithms described in this 
paper have been implemented in Macaulay 2 

1.1. Notation. Throughout we shall denote the ring of polynomials if [xi, .. . ,Xn] 
by if[a;], the ring of polynomials K[di, . . . ,9„] by K[d], and the ring K[x]{d) of 
JsT-linear differential operators on K[x] by D. 

Let us also explain the notation we will use to write maps of left or right D- 
modules. As usual, maps between finitely generated modules will be represented 
by matrices, but some attention has to be given to the order in which elements are 
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multiplied due to the noncommutativity of D. Let us denote the identity matrix of 
size r by id^, and similarly the identity map on the module M by idjv/- 

Given an r x s matrix A — [aij] with entries in D, we get a map of free left 
D-modules, 

where and are regarded as modules of row vectors, and the map is matrix 

•A 

multiplication. Under this convention, the composition of maps — > and 
D" ^ D* is the map — > D* where AB is usual matrix multiplication. 

In general, suppose M and N are left I?- modules with presentations /Mq 
and D'^/Nq. Then the matrix A induces a left D- module map between AI and 

N, denoted [D"^ /Mq) ^ (D^/No), precisely when L ■ A e No for aU row vectors 
L G Mq. This condition need only be checked for a generating set of Mq. Conversely, 
any map of left Z?-modules between M and N can be represented by some matrix 
A in the manner above. 

Now let us discuss maps of right _D-modules. The r x s matrix A also defines a 
map of right D-modules in the opposite direction, 

(Or^iD'-f : [e[,...,Q^^A-[i[,...,Q^, 

where the superscript-T means to regard the free modules {D^Y' ^^^d [D^Y' 
consisting of column vectors. This map is equivalent to the map obtained by ap- 
plying Hom£i(— to D"^ D^, thus [D'^Y' '^^Y alternatively be regarded as 
the dual module HomD(D'',D). We will suppress the superscript-T when the con- 
text is clear. As before, the matrix A induces a right D-module map between right 
D-modules TV' = (D'Y/NI^ and M' = {D^'Y /M^ whenever A ■ L e for aU 

column vectors L £ 7V^. We denote the map by {D'Y/K ^ {D^'Y l^i'o- 

1.2. Left-right correspondence. The category of left D-modules is equivalent 
to the category of right _D-modules, and for convenience, we will sometimes prefer 
to work in one category rather than the other - for instance, we will phrase all 
algorithms in terms of left D-modules. In the Weyl algebra, the correspondence is 
given by the algebra involution 

D^D : x-^dl^ ^ {-dfx°'. 

The map r is called the standard transposition or adjoint operator. Given a left 
_D-module D"^ /M^, the corresponding right ZJ-module is 

^{w)--=^y -(^^o) = Mx)|XeMo}, 

Similarly, given a homomorphism of left D-modules : / Mq — >D'^ / defined 
by right multiplication by the r x s matrix A = [a^], the corresponding homo- 
morphism of right D-modules t(0) : D^/t(Mo) — /t{Nq) is defined by right 
multiplication by the s x r matrix t{A) := [T(aij)Y ■ The map r is used similarly 
to go from right to left D-modules. For more details, sec p^ . 

2. Basic Isomorphism 

The following identification, which we take from Bjork is our main theoretical 
tool to explicitly compute homomorphisms of holonomic D-modules. 
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Theorem 2.1. Let M and N be holonomic left D-modules. Then 
(3) Ext'c (M, N) = Tor,^_, (Ext^ (A/, D),N). 

Proof. Since it will be useful to us later, we give the main steps of the proof here. 
The interesting bit of the construction is the transformation of a Hom into a tensor 
product. The presentation is adapted from Let X* be a free resolution of M, 



: 



•Mo 



M ^0 



We may assume it is of finite length by virtue of Hilbert's syzygy theorem - namely, 
Schreyer's proof and method carries over to D (see e.g. Q). The dual of X* is the 
complex of right ZJ-modules, 



}iomD{X',D) : ^ (D''-)^ 



iceicc a 



degree 



Since YIou\d{D\ D) ®d N Home iV) , we see that nonioiX' ,D)®d N ~ 
HoTao{X' ,N), whose cohomology groups are by definition Ext^(M, iV). Now as 
is customary, replace iV by a free resolution Y' , which we may also take to be of 
finite length. 



(4) y : ^ D'-" ^^^^^ > D'-' ^ D''" ^ TV ^ 

We get the double complex Hom£i(X', D) 0o Y*, 

(5) 



_i -No 





T 

0^ (D''-»)^(giDD=o 



(M_a + l-)»idso 



(Mo 



0^ (D"-")" ®D-D° 



(- idr_^)"®( No) 

(A-f_„+l-)»id3_i 



id,._^ »(-No) 



idrg »(-No) 



(Mo-)®i<is_i 



0^ (D*^ 



(M^„+1-)S 



(Mo-)»ids_(, 



Since the columns of the double complex are exact except for at positions in the 
top row, it follows that the cohomology of the total complex equals the cohomology 
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of the complex induced on the table of Ei terms (vertical cohomologies) , 
(6) 

- Romn{D^-^,N) • ■ • Hom,, (i?- , iV) ^ 

^ v ' " v ' 

degree a degree 

As stated earlier, these cohomology groups are Ext^(A/, N). 

On the other hand, since M is holonomic, the complex Hom£)(X*,Z3) is exact 
except in degree n, where its cohomology is by definition Ext^(M, D). Hence the 
rows of the double complex are also exact except at positions in the n-th column, 
i.e. the column containing terms (£)''-" ®d {—)). It follows that the cohomology of 
the total complex also equals the cohomology of the complex induced on the other 
table of El terms (horizontal cohomologies), which in this case is 

idExt'A (M,D) ®(-JVo) 

(7) ^ ExtJ(M, D)(E)D D"-" ► Ext^(M, D) (gjo ^ 

By definition, the above complex has cohomology groups Tor^(Ext^(Af, D), iV), 
which establishes the identification. □ 



In the next few sections, our goal will be to compute an explicit basis of coho- 
mology classes of the complex (||). In particular, the cohomology in degree cor- 
responds explicitly to Hom£i(M, N) because any map ip G Hom£)(L)'''', N) which is 
in the degree kernel, i.e. in 

(8) H'{RomD{D'--\N) Hom^(i^-° , TV) ^ 0), 

degree 

factors through M ~ 0"^° /A4q, hence defines a homoniorphism ip : M N. The 
reason why it is hard to compute these cohomology classes is that the modules 
liomi:){D^' , N) in the complex (|^) are left _D-modules while the maps Hom£)((Mi-), N) 
are not maps of left Z?-modules. In the next few sections, we will explain how the 
ingredients of the proof of Theorem can be combined with the restriction algo- 
rithm to compute the desired cohomology classes. 



3. Polynomial solutions 



In this section, we give an algorithm to compute Hotlid^M, ^^[x]) for holonomic 
M. This vector space is more efficiently computed by Grobner deformations as 
described in , but we wish to discuss this special case in order to introduce the 
general methodology. 

For N — K[x\, the isomorphism ^ of Theorem |2.l| specializes to 

(9) 



Ext^(M, K[x]) ~ Tor^_,(Ext'^(A/, D),K[x]). 



In this case, the proof of Theorem also leads directly to an algorithm. As a 
U-module, the polynomial ring has the presentation K[x] ~ D/D- {di, . . . , and 
can be resolved by the Koszul complex. 



IC' : 



D 



■ [(-i)"-io„,-,ai] 



D 



degree n 



degree 
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The complex whose cohomology computes Tor^_j(Ext^(M, D),K[x\) then spe- 
ciahzes to Ext2)(M, D) On /C* and is equivalently the derived integration complex 
of Ext2)(M, D) in the category of right D-modules. Oaku and Takayama's integra- 
tion algorithm can now be applied to obtain a basis of explicit cohomology classes 
in iJ"(Ext2,(M, D) ®d IC') ~ Tor^ (Ext^, (M, D), K[x]). These classes can then be 
transferred via the double complex (|^) to cohomology classes in the complex (||), 
where they represent homomorphisms in Hom£)(M, i4r[x]). The method and details 
are probably best illustrated through an example. 

Example 3.1. Consider the Gelfand-Kapranov-Zelevinsky hypergeometric system 
AIa{P) associated to the matrix A = {1, 2} and parameter vector /3 = {5}, i.e. the 
D-module associated to the equations, 



u = 9i + 202 

Here, 9i stands for the operator Xidi. 
A resolution for M^(/3) is 



dl-d2 



while a resolution for K[xi, X2] is the Koszul complex, 

r d2 1 

IC':0~.D i^i^ D ^ 



The augmented double complex Hom£)(X', I?) (Sid IC' is 



ExtS(M^(/3),i?) 



K[xi,X2] 



K[xi,X2] 



l-v u+2]- 



d2 1 



d2 1 



92 
82 

-di 
-di 



d2] 



Ext^(MA(/3),i?) 



-V u+2 




1 

— u u+2\ 



■[di 82] 



[ 9i 82 



[-D U + 2]- 



K[xi,X2] 



[-%] 



■[di 82] 



Here, we interpret an element of a module in the above diagram as a column 
vector for purposes of the horizontal maps and as a row vector for purposes of the 
vertical maps. The induced complex at the left-hand wall is the derived integration 
to the origin of Ext_D(MA(/3), D) in the category of right D-modules. Applying 

the integration algorithm, we find that the cohomology at the module 



in the 



bottom left-hand corner is 1-dimensional and spanned by the residue class of 



^1,0 



-{2x1x2 



AOxlxl 



120x1X^)51 



(4 



30x1x2 + 180x^X2 - 120x-^ 
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We lift this class to a cohomology class of the complex induced at the top row via 
a "transfer" sequence in the total complex given schematically by 



u+2 








d2 " 
d2 
-di 
-di 



1,1 



■[9i ^2 



9 Li,o 

In other words, Li i is obtained by taking the image of Li.o under the vertical map 
and then a pre-image under the horizontal map, and similarly for Li^2- We find 
that, 



Li. 



2x\x2 
-{x\ - 2Qxlx2 
- {xl - 2{)x\x2 
{x\ - 20x1x2 4 



+ I20xixl 
f 60xix|) 

f mxlxl) 

60x1X2)91 - 



(lOxf - l2Qx\x2 



l2Qxl) 



[ 



20xiX2 + GOxix^ 



The space of polynomial solutions is spanned by the residue class of Li^2 in K[xi, X2], 
which is xf — 20xiX2 + GOxiXj. 



Remark 3.2. The transfer sequence above is used to show that Tor is a balanced 
functor in Weibel (l8| . A generalization of the transfer sequence is also used by the 
second author to compute the cup product structure for de Rham cohomology of 
the complement of an afline variety in . 

From a practical standpoint, the method outlined above is not quite the final 
story. The detail we have left out is how Oaku and Takayama's integration al- 
gorithm actually computes the cohomology classes of a Koszul complex such as 
Ext2)(A/, Z?) JC*. Their algorithm does not compute these classes directly. 
Rather, their method (phrased in terms of right D-modules) is to first compute 
a V-strict resolution Z' of Ext^(Af, I?). Then they give a technique to compute 
explicitly the cohomology of Z* iS)d K[x]. This complex is quasi-isomorphic to 
Ext^(M, D) K*, and cohomology classes are transferred by setting up another 
double complex Z* IC* . Thus, our method as described to compute polynomial 
solutions requires two transfers via two double complexes. 

Given the true nature of the integration algorithm, the two transfers can be 
collapsed into a single step. Namely, we start with Hom£)(X*, 13), 

Homz3(X',i3) :0^.-. CP^-^f . . . ^JH^ (Ji^ ^ Q 

degree n degree 

which is exact except in cohomological degree n because M is holonomic. We are 
interested in explicit cohomology classes for iJ°(Hom£)(X*, D)(g)£)/C[a;]). To obtain 
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them, we replace }iom£){X* , D) with a quasi-isomorphic l^-adapted resolution E* 
along with an explicit quasi- isomorphism tt, from E* to HomoiX' , D). That is, 
we make a map 7r„ from a free module (D'*"")^ onto some choice of generators 
of ker(M_„-), take the pre-image P of im(M_„+i-) under Wn, and compute a V- 
adapted resolution E' of _D*""/P. Schematically, 



0^...^(_D'- 



Mo- 



Using the integration algorithm, the cohomology classes of the top row can now be 
computed. In order to transfer them to Ilomj){X* , D) (g)^) IC[x], a chain map lifting 
7r„ is computed and utilized as suggested by the dashed arrows. We summarize the 
algorithm as follows. To keep computations in terms of left -D-modules, we make 
use of the transposition t at various places. Applying r to the polynomial ring 
gives the top differential forms H. — D/{di, . . . , dn} — t{K[x]). 

Algorithm 3.3. [Polynomial solutions by duality] 

Input: {Li, . . . , LrJ C D''" such that M = {D^'o /D ■ {Li, . . . , L^}) is holonomic. 
Output: The polynomial solutions R G if [x]""" of the system of differential equa- 
tions given by Li • R — 0, i ^ I, .. . , ri. 

1. Compute a free resolution A"' of M of length n + 1. Let its part in cohomo- 
logical degree —n be denoted: 

•M_„ -M- 



Form the complex r(IIomc(X', £))) obtained by dualizing X' and then ap- 
plying the standard transposition. Its part in cohomological degree n now 
looks like: 



•r(A/_„) 



Compute a surjection 7r„ : D''-" ker(-r(M_„)), and find the pre-image 
r(P) := 7r„-i(im(-T(M_„+i))). This yields the presentation D^~"/t{P) ~ 
r(ExtS(M,i^)). 

Compute the derived integration module (i^j^ {D'^-" /T{P))[n]) using 

Algorithm 8.S. In particular, this algorithm produces 

(i.) A F-strict free resolution E' of D^-"- /t{P) of length n + 1, 

£" : ^ ^ 0"-"+^ < ^ D"-^ ^ ^ D'K 

degree n degree 

(ii.) Elements {gi, . . . ,gk} C -D^" whose images modulo im{il D''^) form a 



basis for 



H' 



E-) 



ker(t70r>D°-'-^O(8j£)°o) 



5. Lift the map 7r„ to a chain map tt, : E' 
maps TTi : D"-^ — > D^-'. 

6. Evaluate {T(7ro(5i)), . . . , T(7ro(5fc))} and let {i?i(x), 
ages in (D/D • . . . , ~ K[^Y°. 



T(IIom£i(X*, Z?)). Denote these 
,i?fc(x)} be their im- 
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7. Return {i?,i(x), . . . , i?fc(x)}, a basis for the polynomial solutions to M . 

Example 3.4. Let us return to the GKZ example and apply the revised algorithm. 
For Step 1, we have already computed a resolution X* . Its length equals the global 
homological dimension. Thus for Step 2, we get a complex which is a resolution for 
the holonomic dual, 

T(Hom(X-, D)):Q^ ^ 0, 

where u' = —{Oi+ 202 + 6) and v' — ~{df +92). For Step 3, it follows that we get 
the presentation, 

riE.tn{M,D)) ^ ^ D ■ {9, + 29^ + 6, + d,} 

For Step 4, we compute the derived integration of this module. It turns out the 
complex r(Hom(X*, D)) is already a V^-strict resolution when taken with the shifts 
^ ^^[O] <- £'^[1,0] ^ D'^[l] ^ 0. The integration 6-function is s - 4, hence 
according to the integration algorithm, (E)d T(llom{X' , D)) is quasi-isomorphic 
to the finite-dimensional subcomplex, 

[F^ and F^ are explained in the appendix.) Here, i^'*(r2^ [— 1]) is spanned by the 
21 monomials of degree < 5, 

{1, Xi, X2, . . . , X]^, X-^X2^ '^1*^2' ^l'^27 -^1*^2' ^2}' 

while F'*(r2^[0, —1]) is spanned by the 36 monomials 

{1, Xi, X2, ■ ■ ■ , Xi, X1X2, X1X2, X1X2, X2} " Gl U 

{1, Xl, X2, . . . , X-^^ X^X2^ ^1^2^ ^1^2^ '^l'^2' '^2} ' ^2- 



The matrix 



induces a map between them whose kernel is spanned by the 



degree 5 polynomial Ri = (xf — 20x^X2 + GOxix^)- 

4. Rational Solutions 

A duality algorithm to compute the dimensions of Ext^(M, if[x][/^-'^]) for holo- 
nomic M was given in ||l^ . In this section, we show how to extend this algorithm so 
as to compute an explicit basis of HomD(Af, ^[x] [/~^]). The method is essentially 
the same as the algorithm for polynomial solutions. Also, since any rational func- 
tion solution has its poles inside the singular locus of Af , we obtain an algorithm to 
compute the rational solutions of M . Finally, we remark that a different algorithm 
to compute rational solutions based upon Grobner deformations was given in Jl2| . 
Here, as otherwise, we shall use N[f~^] to denote N C15k[x] K[x][f^^]. 

Algorithm 4.1. (Rational solutions by duality) 

Input: {Li, . . . ,LrJ C D''° such that M = {D''°/D ■ {Li, . . ■,Lr,}) is holonomic. 
Output: The rational solutions R G (x)''° of the system of differential equations 
given by Li • R — 0, i — 1, .. . , ri . 

1. Compute a polynomial / which defines the codimension 1 component of the 
singular locus of M (see e.g. Jl5|). 
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2. Compute a free resolution X* of M up to length n + 1. Let its part in 
cohomological degree —n be denoted: 

. . . _> O^-^-^ £)'■-'« '^'^~" + ; D^-n+l _5, . . . 

3. Form the complex T(Hom£)(X*, 13)). Its part in cohomological degree n now 
looks like: 

4. Compute a surjection 

Wn : ^ ker(-T(Af_„)), 

and find the preimage t{P) :— t37„^^(im(-T(Af_„+i))). Denote by tn^„ the 
induced map on the localizations, [/^i] ker(-T(M_„))[/-i]. 

5. Compute the localization of D'^-" /t{P) at / using the algorithm of ||T^. This 
produces a presentation. 



-iQ) \r{P)[f-' 

a mod t{Q) — ^ (ci mod r(P)) ® /~°* 
6. Compute the derived integration module H^{{^1 D*-"/r(Q))[n]) using 



Algorithm 8.S. In particular, this algorithm produces 

(i.) A y-strict free resolution E* of D'^-" /t{Q) of length ?i + 1, 

E' : ^ ^ 0^"-"+^ < ^D"-!^ ^ D'K 

degree n degree 

(ii.) Elements {51, . . . , gk} C whose images form a basis for 

7. Let 7r„ be the composition 

tx7;,„ o ^ : D'- ^ I3^-"[/-i] -> ker(.T(M_„))[/"'], 

where (p : D^-" — > [f~^] is the map defined by 1-^ /""'. Lift 7r„ 
to a chain map tt, : _B* — > IIomc(-f , -D)[/"^]- 

8. Evaluate {r(^o(5i)), • ■ • , C i3'^"[/-i] and let {i?i(x), . . . , i?fe(x)} 
be their images in iD[f-^]/ D[f-^] ■ {Oi, . . . ,a„})'-« ~ A'[x][/-i]'-°. 

9. Return {i?i(x), . . . , i?fe(x)}, a basis for the rational solutions to M. 

Proof. As explained in , any rational solution of M has its poles contained inside 
the singular locus of M. The proof is now essentially the same as for the polynomial 
case. The space of rational solutions can be identified with the 0-th cohomology of 
the complex (^, which specializes to 

~ RoinD{X',D)[f-'] D[f-'] ®D K[x] 

~ EomDiX',D)[f-^](^DK[x] 

Since the complex HoiaiDiX' , D) is exact except in cohomological degree n 
where its cohomology is Ext^ (M,D), and since localization is exact, the com- 
plex Homo {X *, D)[f~^] remains exact except in cohomological degree n where 
its cohomology becomes Ext2)(M, D) [/"!]. Hence BomoiX' , D)[f-^] i^id K[x] 
computes the derived integration modules of Ext^{M, D)[f^^] in the category 
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of right Z?-niodules. The above algorithm computes cohomology classes for the 
derived integration modules and transfers them back to cohomology classes of 
HomD(X',D)»z5if[x][/-i]. □ 

Remark 4.2. Let us explain how the lifting of 7r„ to a chain map in Step 7 
may be accomplished algorithmically. We wish to do computations in terms of 
D and not D[f^^]. The idea is that localization is exact, hence any boundary in 
Hom£i(X*, D)[f^^] is the localization of a boundary in HouiuiX* , D). Suppose we 



have computed tt, : 



Then to compute T^j-i, we first compute 
the images £i of G D^j+i under ZJ'^-i+i D''-' — ^ D^-i[f^^]. Because the 
existing 7r„, . . . ,7rj are the beginning of a chain map, the ii are in the image of 
D'^-^+^[f~^] — > D'^-^[f~^\. Now we use the fact that localization is exact, which 
means for sufficiently large mi, f™^£i is in the image of D^-^+'^ D^^^ . To find 
valid rrii, we can multiply £i by successively higher powers of / and test for mem- 
bership at each step via Grobner basis over D. Now compute any preimage Pi of 
in D^-^ . The map nj^i : — > D^^^^^[f^^] may be defined by sending 

6i ' ^ ym; Pi- 



Example 4.3. The following system of differential equations of two variables is 
called the Appell differential equation Fi{a, b, b',c): 

OA^x + 9y+c-l)- x{e, + Oy + a){e, + 6), 

ey{e, + ey + c~i)^ y{0, + 9y + a){ey + b'), 

{x - y)dxdy - b'dx + bdy, 

where a,b,b',c are complex parameters. In [|l2|, the dimension of the rational 
solution space of i^i(2, — 3, — 2, 5) was computed using the duahty method. This 
system has rank 3, and its solution space is spanned by a polynomial, a rational 
solution with pole along x, and a rational solution with pole along y. 

Let us obtain the solution with pole along x explicitly. In [ p2[ was computed a 
resolution for Fi{2, —3, —2, 5), 



0, 



so that T(Hom£,(X',D)) is a resolution for T(Ext^(Af, D)) = 4, 2, -3), 

0. 



r{Romn{X',D)) : Q^D' 



where, 
r(M-i) 

T(Afo) 



{9., + \)dy - {0y + 2)9, 
(y2 _ y){d^dy + dl) + 2{x + y)d^ - 2ydy ~ 29, + 79^-2 
(y' - y){d.dy + dl) + 2{x + 2y)d, ~ 39, + 69^ 



V 



-(0,+4)9y + (&y + 3)9. 
It was also computed that the localization has presentation 
D ~ / D[x-^] 



r{Q) 



im{-T{M-i))[x- 



(ExtS(Af,i?))[x- 



1 mod t{Q) 



(1 mod im(-T(M_i))) ® a;"'^ 
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where 



t{Q) = D 



(O^Oy +01+ 89y + 26, + 12) - (6, +0y+ A)dy 



{OJy + 20^ + 70y + 14) - + 10)xdy 

and that D/t{Q) has a F-strict resolution, 



E' :0 



where 



= x^dxdy — xydxdy — 2xdx + llxdy — 7ydy — 14 
U2 = x^dl + xydxdy - x^dl - x^d^dy + IQx^dx + nx^dy+ 
4:xydy — 9xdx — llxdy + 52x — 7 
= x^d^ + xydxdy - x'^dl - xydxdy + IGx'^dx + 12x'^dy+ 

Axydy — Sxdx — llxdy + h2x — 6 
= —x'^dxdy + xydxdy + 2xdx — llxdy + 6ydy + 12. 



Ml 



Vi 



V2 



We would like to construct a chain map tt, : £* — > T{liomn{X* , D)) which 



lifts the map ttj : D^[0] D[x-^]'^ defined by 1 



To compute the 



next map tti : -D^[0, — 1] D[x ^]^, we need to find preimages of the elements 
772 ° (•['"1, W2]^)(ei)and n2 o (•['Ui,'U2]'^)(e2)under the map (•t(M_i)) : Z)^[a;~^] < — 
£)2[a;-i]. Note that 

7r2 o (•[ui,U2]'^)(ei) = ui-a;-^ 

= x-%{ex + 4.)dy-iey + 2)dx) 

It follows that 7r2 o {■[ui,U2]^){x^ei) = (•r(M_i))(ei) so that we may set 7ri(ei) = 
x~^e'i. In a similar manner, we obtain the chain map. 



E* : 



T(Homr,(X*,£))) : 

where 

a 
b 

c 



D^[x-^] 



D'[0,-l] 



D^[x-^] 



■IVL,V2] 



D'[-l] 







•t(Mo) 



D^[x-^] ^0, 



-^{y^ -y){dx + dy)+x + 2y 
^{x-y)dx + 2 
l{x-y)dx + ^ 



The integration 6-function is (s — ll)(s— 4)(s — 1), hence according to the integration 
algorithm, fl E* is quasi-isomorphic to its subcomplex F^^{Q. 0£) E*). Using 



Macaulay 2, we find that ker(f] 
fli 



f2^) is 2-dimensional and spanned by, 



x»?y - |.T» - 6.T«y + f .t8 + ^x''y 



_ 126^7 _ M^6y + i|8^6 + ^^5y _ 21x5 
92 = -x^ + f x^y + ^x"^ - f^xy - + ^y 
The residue class of T(7ro((7i)) yields the polynomial solution, 

, „ 2 144 165, 2 
(-9y2 + —y - _)x2 



24 

(2y2-6y+— )a;-'^ 



,72 2 252 252, , 42 ^ 168 
+ {^y^--^y+-^)x + {--y-' + —y-A2), 
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while the residue class of t (710(32)) yields the rational solution, 

/i4 ^3 ^22 ^3 2 ^2^2^ ^2 

-6a; +4xy--xy + 4x ~ —x y + -xy - -x + -xy ~ —y 
By similar methods, we obtain the rational solutions with pole along y, 

^3^2 _ ^^2^3 + 21 _ ,,y, _ 5^3y + Q^2y2 X _^ 
-^^Xy^ + ly^ + ^^X^-^lx^y + ixy^-ly^ ' 

Together these solutions span the holomorphic solution space in a neighborhood of 
any point away from xy = 0. 

Remark 4.4. In the next section, we give an algorithm to compute Hom/)(Af, N) 
for arbitrary holonomic M and N. Using it with N = _ft'[x][/^^], we get a similar 
but computationally different duality method to compute rational solutions. The 
basic difference is that the algorithm of this section uses computations over D and 
in principle over D[f~^], while the algorithm of the next section uses computations 
over the Weyl algebra in twice as many variables. From the computational 

perspective, we believe the algorithm of this section is more efficient. 

5. Holonomic solutions 

In this section, we give an algorithm to compute a basis of Hom/j (M, N) for 
holonomic left U-modules M and N . We will use the following notation. As 
before, D will denote the ring of differential operators in the variables xi, . . . ,a;„ 
with derivations 9i, . . . , i9„. Occasionally we will also write or for Z3. In a 
similar fashion, Dy will stand for the ring of differential operators in the variables 
yi, ■ ■ ■ ,yn with derivations . . . , i5„. 

If X is a Da^-module and Y a D^-modulc then we denote hy XMY the external 
product of X and Y . It equals the tensor product of X and Y over the field K , 
equipped with its natural structure as a module over = Kl Dy, the ring of 
differential operators in xi, . . . , x„, j/i, . . . , j/n with derivations {di, i5j}i<ij<n- In 
addition, let rj denote the algebra isomorphism, 

. ^ ^ { Xi ^ \xi - S^, di (-> ij/j + di 

V- 2n * 2n \ y, ^ -^Xi - 5^, 5^ ^ \y^ - J .^^ 

and let A and A denote the right _D2n-inodules, 

A:=7 ^ ^ ^ A:= ^ = r;(A). 

{Xi - yi,Oi + di:l<l<n} ■ D2n ^D2n + yD2n 

As mentioned in the introduction, an algorithm to compute the dimensions of 
Ext^ (M, A'') was given in based upon the isomorphisms (|l|) and (||) : 

Ext^(M, A^) = Tor^_,,(ExtS(M, D),N) 

Torf (M', A^) ~ ToTf^-{D2n/{x^ ~ y„ 9. + S,}^^, ■ i?2«, r(M') K N). 
Combining these isomorphisms where M' = Ext^) {M, D) produces 
(10) 

Extl5(Af, N) ~ Torf" {D2n/{x^ -y^,^, + S,}^^, ■ D2,., r(Ext£(M, D)) M N) 

In order to compute HoiiidIM, N) explicitly, we will trace the isomorphism (|lO|). 
We explain how to do this step by step in the following algorithm. The motivation 
behind the algorithm is discussed in the proof. 
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Algorithm 5.1. (Holonomic solutions by duality) 

Input: Presentations AI = /Mq and N = D^" /Nq of holonomic left D-modules. 
Output: A basis for HomD(M, TV). 

1. Compute finite free resolutions X* and Y' of M and N, 

■M-a + l 



A* : ^ 

degree —a 



Y' -.0 



D 



degree 



M 



A ^ 



degree —h degree 

Also, dualize A* and apply the standard transposition to obtain, 



r(Homc(A*, D)) : ^ 

degree a 



degree 



2. Form the double complex T(Homo(A*, D)) Kl Y* of left r'2n-modules and its 
total complex 

Z' -.0^ D2n^ < ^ D2r^ < ^ D2n'-' ^ 

degree a degree 

where 

i-j=k 

Let the part of Z* in cohomological degree n be denoted, 



•T„_i 



D2n 



3. Compute a surjection 7r„ : £'2n""-» ker(-77(T„)), and find the preimage P := 
TT-i (im(.r?(T„_i))). 

4. Compute the derived restriction module i/"((A®£^ -D2„""/P)[n]) using Al- 
gorithm ^.6| . In particular, this algorithm produces, 

(i.) A y-strict free resolution E* of D'^^ / P of length n + 1, 

E':Q^ D2n^ ^ i?2„""-^ ^ . . . ^ 7^2„"^ - D2n^ ^ D2„"-^ 
degree n degree 

(ii.) Elements {gi, . . . , g/c} C D2n^° whose images in A ®d2„ E* form a basis 
for 

770 ((a ®^ Ii^)[„]Ui70(A^^^ j^.)^ ker(A«.,„D.„-i^A^„„.I.. "0) 

5. Lift the map 7r„ to a chain map tt, : —>■ rj{Z*). Denote these maps by 

6. Compute the image of each gi under the composition of chain maps. 



A^D^^Z' Tot*(HomD(A*,D) (g)DY*) 



pi 



viz') 



Hom£,(A', A) 
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Here pi is the projection onto HoniD(X*,£') eg) y° followed by factoriza- 
tion through Nq. These are all chain maps of complexes of vector spaces. 
Step by step, we do the following. Evaluate {Li = r]~-^{TTo{gi)), . . . , Lk = 
r]~^{'^o{gk))}, and write each Li in terms of the decomposition, 

j 

Now re-express Li^ modulo {xi — yi,di+6i : 1 < i < n} ■ D2n®D2„{D^°^D'^°) 
so that Xi and dj do not appear in any component. Using the identification 
jjro ^ jjso ^ ei © ■ • • ® DanSr, whcrc {e,} forms the canonical D-hasis for 
^ we then get an expression 

Let . • . ,?i,ro} be the images in {D^°/Nq) ~ A^. Finally, set G 

Hom£)(Af, A^) to be the map induced by 

{ei 4,1, 62 ^ ^i,2, . . . , Cro 1-^ ^i,ro}- 

7. Return . . . , a basis for Hom£)(Af, N). 

Proof. The main idea behind the algorithm is to adapt the proof of Theorem 2.1, 
In that proof, we saw that Tot'(Hom(X',i:))(8)Dy) -^IlomD{X',N) IS a quasi- 
isomorphism. Thus it suffices to compute explicit generating classes for 

H°{Toi'{Ylo\T^D{X' ,D) ®Dy')) H"{}iomD{X',N)) ~ UomoiM, N). 

Here, the double complex Hom£)(X', D) (g)^ Y* is in some sense easier to digest 
because it consists entirely of free ZJ-modules. However, it too only carries the 
structure of a complex of infinite-dimensional vector spaces, making its cohomology 
no easier to compute than the cohomology of KomoiX* , N). 

Thus, we instead are led to consider the double complex r(Hom£)(X*, D)) Kl Y* 
of Step 2, whose total complex T* does carry the structure of a complex of left D2n- 
modules. Moreover, we can get back to the original double complex by "restricting 
back to the diagonal" . In other words, we claim that as a double complex of vector 
spaces, Hom£)(X*, D) (E)dY' can be naturally identified with the double complex, 

A(g,D {TiRomDiX',D))MY'). 

To make the identification, first note that the natural map 



^ {xi - yi,di + 5i ■■ I < i < n} ■ D2n 

is an isomorphism of left I?j,-modules. Let {ei, . . . , e^} denote the canonical basis 
of a free module . Then an arbitrary element of A ®D2n i^x^ ^ ^y^) '^^^ be 
expressed uniquely as X]fe ^fc ^ '^fe' where e Dy'' . Similarly, an element of 
iSidD^ can be expressed uniquely as ek^mu where ruk e D'^ . Hence we get 
an isomorphic identification as D^-modules of /^®02n {Dx^ ^ Dy^) and ®o D" ■ 
In particular, this shows that the modules appearing in the double complexes are 
the same. 

It remains to show that the maps in the double complexes can also be identified. 
An arbitrary vertical map of A ®D2„ (T(HomD(^*, D)) M Y*) acts on an arbitrary 
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element J2k 1 ® efe ^ f^k according to, 



idA<8i(-id^i)*KI(-JV^) 



This is exactly the way the corresponding vertical map in HomD(X',D) Y* 
works on the corresponding element: 



(-id.,)'®(-Arj) 



Likewise, an arbitrary horizontal map of A <8>d2„ (T(Hom£)(X*, £))) IE Y*) acts 
on an arbitrary element according to, 



Here, we would like to re-express the image 1 (8) (■T(Mi))(efc) Kl mfc in the form 
J2k If^CfelElnfe. To help us, note the following computation in A^/j^^ (Dx^^Dy^): 

(1 (8) aj^a'^ei Kl m) = 1 O a'^Cj Kl y"m = 1 {-6fy"m = 1 K T{y"'5^)m. 

Using it, we get that 

k k j 

= ^^lOCj KT(r(Mi)jfe)mfc 

k j 

= ^ ^ 1 e^- Kl {Mi)jkmk 

k j 

This is exactly the way the corresponding horizontal map in Hom£)(X*, D) ^dY* 

works on an arbitrary element: 

iMi-)<S,ids, 

L>'^-+i «n D'^ '—^ £)'■• 0n D'^ 



Efcefcigimfc Efc Ej «3®(-'^»)jfe™fc 

Thus, we have given an explicit identification of A^p {T(Romjj{X* , D))^Y*) and 
RomD{X\D) (^dY*. 

The task now becomes to compute explicit cohomology classes which are a basis 
for H^{A <8)D2,i Z*). To do this, we note that Z* is exact except in cohomologi- 
cal degree n, where its cohomology is r(Ext2)(M, D)) M N. This follows because 
T{iloTaD{X* , D)) is exact by holonomicity except in degree n, where its cohomol- 
ogy is T(Ext^(M, D)), and Y* is exact except in degree 0, where its cohomology is 
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N. In other words, the complex A ^Ua^ Z* is in some sense a restriction complex. 
Namely, after applying the algebra isomorphism rj, we get an honest restriction 
complex A (g) rj{Z') for the restriction of ?7(r(Ext2)(M, D)) M N) to the origin (the 
restriction complex of a left I?2n-module M' is by definition A <E)i)^^ M'). 

We can thus compute the cohomology groups of A ®d^„ v{Z*) by applying the 
restriction algorithm. However, since we are after explicit representatives for the 
cohomology classes, we need to use a presentation of ?7(r(Ext^(M, D)) MN) which 
is compatible with r]{Z*). This is the content of Step 3. Once equipped with 
a compatible presentation, we apply the restriction algorithm to it, which is the 
content of Step 4. This step produces explicit cohomology classes of A C5i£)2„ E' , 
where E* is a F-strict resolution of r]{T{Ext1){M, D)) Kl N). To then get explicit 
cohomology classes of A ®d2„ v{Z*), we construct a chain map between E* and 
ri{Z*), which is the content of Step 5. The cohomology classes can now be trans- 
ported to A (8)D2„ viz') using the chain map, then to A ®£)2„ Z' using ?]~^ , then 
to Tot*(Hom£)(A'*, D) (giu Y*) using the natural identification described earlier, 
and finally to the complex Hom£i(A"', A^) using the natural augmentation map. 
These steps are all grouped together in Step 6. This completes the proof of the 
algorithm. □ 



Remark 5.2. In Algorithms 3.3 and 4.1 we used the integration algorithm as the 



main workhorse, while in 5.1 we used restriction. As should become apparent from 
the appendix, these are really mirror images of each other. Rational and polynomial 
solutions naturally fit into the integration picture. On the other hand, most papers 
are written in the language of restriction. 



Example 5.3. Let M = D/D ■ {d - I) and N = D/D ■ {d - 1)^, where D is the 
first Weyl algebra. Then for Step 1, we have the resolutions. 



X' : 



^liO^^i 







r* : 







For Step 2, we form the complex Z* = Tot(r(Homz5(X', D)) M V) 

■ [(8„-l)^-(a,+l)] 



Z' : 



D2' 



degree 1 

For Steps 3-5, we get the output. 



D2^ 

dcercc —1 



viz') : 0^ D2' 



(^y-Ox-i)- 



Do 



77l = -[ll 



E' : 



degree 

'^o=-||!y-k-i 1] 



^0 



-1,2] 



•[y^-^y-Ox-l] 



The complex E* is a V^-strict resolution of the cohomology of at degree 1, and 

the restriction &- function is = (s + l)(s + 2). Hence A(SidE' is quasi-isomorphic 
to its sub-complex F~-^(A ®d S*) 

'^y+Ox+il r 0©T 

0®^ 



0^0 



Span 



K 



Span^{l} ^ 



o©a, 
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Hence the cohomology H^{A (S)d E') is spanned by {0 © 1,0 © dy}. Applying 
ttq, H^{A (g)£i r]{Z')) is spanned by the images of — dx — I) (B l,dy{^y — 
dx — ^) ® 9y}. Next applying r]~^, i?°(A Z') is spanned by the images of 
{Li = {dx + 2dy-l)®l,L2 = -^{xdx + 2ydy + ydx + 2xdy~x~y)®~^{x + y)}. 
Modulo the right ideal generated by {x — y,dx + dy}, we can re-express these 
cohomology classes by {{dy — 1) © 1, {ydy — y — 1) © —y}. Applying pi we get 
{Lifi = dy ~ I, L2fi — ydy — ?/ — 1}, which corresponds to a basis of HomD(M, N) 
given by, 

D .[d-i] D 



D-{d' 
D 



I) 



■[(xd- 



D-{d- 
D 



I) 



D-{d-i) 



D-{d-lY' 



6. Extensions of holonomic £>-modules 

In this section we explain how one can modify our algorithm for the compu- 
tation of Hom£)(Af, A^) in order to compute explicitly the higher derived functors 
Ext^(M, N) for holonomic D-modules M and A^. 

A useful way to represent Ext^(M, A) is as the i-th Yoneda Ext group, which 
consists of equivalence classes of exact sequences, 

-i+2 ^ ^ ^ ^ ^j^O 



?:0 



A^ 



Q 



M 



0, 



for any list of (not necessarily free) 13- modules Q, A . . . Two exact se- 

quences ^ and ^' are considered equivalent when there is a chain map of the form, 

^ : — > A Q A"*+2 — , M — > 




M 



0. 



In our modified algorithm we follow the same steps as in Algorithm 5.1, except 



that in Step 4 we compute iJ^"+*(A (£'2n""/-P)) instead of i/-"(A (g)^^ 



{D2n^" / P)) . The output is a basis {931, . . . , tp/c} of the finite-dimensional AT-vector 
space iJ*(IIom£)(A*, A^)), where X' is a free resolution of M, 

■Mo 



X' -.Q 



degree —a 



■M- 



degree 



M ^ 0. 



To obtain the i-th Yoneda Ext group from our output for Ext^(M, A^), we 
follow the presentation of [^8[ Section 3.4] and associate to a cohomology class 
(f G iJ*(IIom£)(A', A^)) the exact sequence, 

£,{ip) : — >N — >Q — > — > > D''" ^ Af ^ 0. 

Here, Q is the cokernel of (-M-i+i, Lp) : D^-' — > D^-^+^ (B N, and the maps are all 
the natural ones. It is worth pointing out that A^ is indeed a submodule of Q by the 
following argument. G 7J*(Hom£)(A', A^)) is computed as a map from A^' to Y'^ 
with the property that — > A~* ^ ^ A^ is zero. Assume a e ker(A^ Q). 
Then (0, a) G im(-A/_i+i, Lp), so there is & e A^* with h ■ M-i+i = and (f{h) — a. 



Since ker(-M_i+i) = im(-M_ 



M-i and it follows that a — ip{c ■ M- 



0. 
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Notice that the only difference between any ^{ip) and ^{ip') are their correspond- 
ing Q's and the maps to and from them. In terms of the basis {(pi, . ..,(pk} of 
i/*(Hom£)(X*, iV)), the set of possible Q's which appear can be packaged as the 
set, 



(0 © TVo) + • {(-M-, + ELi ^hV>h){ej)V-{' 



(Ki,...,Kfe) e k'' 



When i — 1 for example, the 1-st Yoneda Ext group consists of equivalence 
classes of extensions of Af by A^, 

^:0-^A^ — >Q — > M = {D'"> • Afo) ^ 

where Q = X'^ © N modulo X''^ ■ A//_i. Thus, once we have computed a basis 



{(fii, . . . , ipk} of H^(HomD{X' , N)) via the modified Algorithm 5.1, the possible 
extensions Q of Af by A*" are, 



Example 6.1. Let D = K{x,d) be the first Weyl algebra and M = D/D ■ d, 
N = D/D ■ X. Then for Step 1 of Algorithm 5.1, we have the resolutions, 



X' -.Q^ D^ D^ ^0, y : ^ L'l ZJi ^ 
For Step 2, we form the complex Z* = Tot(r(HomD(X', D)) K y), 

2 -[y^-d.] 



Z' -.0 



degree 1 



D2 



Do 







For Steps 3-6, we find that H^{A Z') is spanned by {1}, and projecting by 
Pi, 'Ext\){M,N) ~ (Y{o\T\]j{X* , D / D ■ x) is spanned by the natural projection 
If : D ^ {D/D ■ x). For k E K, the cohomology classes nip correspond to the 
extensions on the bottom row of the following diagram, 

D 











D 



D 



D 



D-d 







[1,0] 



■[0,1] 



£> • el © D • 



id 



D/D-8 



D 

Dd 







D ■ x D ■ xel + D ■ {kbi + de'2) 

When K 7^ 0, the module Q{k) = {D ■ e{ ® D ■ e^)/(£> • xe{ + D ■ (reel -I- 9e5)) is 
generated by €2 and is always isomorphic to D/D ■ xd. When k = 0, the module is 
no longer generated by 62 and is not isomorphic to D / D ■ xd. 

In fact, the module {D-e{®D-e2)/{D-xei + D-{nei+de2)) is always generated by 
the residue class of el -I- e'2 and has the cyclic presentation D/D ■ {d^x + nxd, x^d} 
with respect to this generator. Using this presentation, the extensions take the 
form, 

D -[-xd] D -[xd+i] D 







D-x 



D ■ {d'^x + Kxd, x'^d} 



D-d 



0. 



20 



HARRISON TSAI, ULI WALTHER 



One can picture Q{k) as the K[x]-mod\ile K[x] + x ^K[x ^] with the twisted 
muhiphcation rule x ■ (x^^) — k which is a direct sum if k = 0. 

7. Isomorphism Classes of £>-modules 

In this section, we give an algorithm to determine if two holonomic _D-modules AI 
and N are isomorphic and if so to produce an explicit isomorphism. For M — N, we 
also give an algorithm to find all isomorphisms from Al to M and mention some well- 
known applications of the endomorphism ring End£)(Af). Here, Endi3(M) denotes 
the space of endomorphisms of a D-module M = D"^^' / Im, where endomorphism 
means Z?-linear maps from M to A/I. Similarly, Isod(M) denotes the units of the 
ring EndD(Af). 

If holonomic AI and N are isomorphic, then HomD(Af, A^) ~ End£)(Af) is a 
finite-dimensional ii'-algebra. In the theory of finite dimensional _fC-algebras, the 
Jacobson radical J is the intersection of all maximal left ideals of E, and it has the 
property that the quotient E/ J is a semi-simple ii'-algebra. By the Wedderburn- 
Artin theorem, a semi-simple algebra is isomorphic to a direct product of matrix 
rings over division algebras, and hence by taking the algebraic closure, we find 
that E/ J ®K K is isomorphic to a direct product of matrix rings over the field K. 
One consequence of this decomposition is that the non-units of E/ J ®k K form a 
determinantal hypersurface. In particular, the units of E/ J ®k K form a Zariski 
open set, and hence the units of E / J also form a Zariski open set. Moreover, units 
and non-units respect the Jacobson radical in the sense that if j is in the Jacobson 
radical of E and if u is a unit of E then u~\- j is also a unit, and similarly, if n is not 
a unit of E then n + j is not a unit. We can thus conclude the following lemma. 

Lemma 7.1. Let M be a holonomic D-module. Then the space of D-linear isomor- 
phisms IsO£)(M) from AI to itself is open in End£)(M) under the Zariski topology. 

The lemma says that if holonomic AI and N are isomorphic then most maps 
from AI to N are isomorphisms. We now give an algorithm to determine if AI and 



N are isomorphic based on Algorithm p.l| and Lemma 7.1. 



Algorithm 7.2. (Is AI isomorphic to A^?) 

Input: presentations A4 ~ D""" /D • {Pi, . . . , Pa} and N ~ D"" /D • {Qi, . . . , Qt} 
of left holonomic D-modules. 

Output: "No" if AI N; and "Yes" together with an isomorphism l : AI ^ N if 
AI ~ N. 

1. Compute bases {si, . . . , Sa} and {ti, . . . , t^} for the vector spaces V = IIom/)(Af 



and W = IIom/)(A^, M) using Algorithm 5.1, where Si and tj are respectively 
rriM X mjv/ and mpf x m^v matrices with entries in D representing homo- 
morphisms by right multiplication. Recall that we view D™'^' and D™" as 
consisting of row vectors. If cr 7^ r, return "No" and exit. 
Introduce new indeterminates {fJ.i}\ and and form the "generic homo- 

morphisms" J^it^i^i ^ IIom£)(Af, iV) and J^j ^j^j ^ HomD(A^, A/). Then the 
compositions J2i j f-i'^j^i '^3 : AI ^ N ^ AI and J2i j t^i'^j^j ' Si : N ^ 
AI —> N are respectively rriM x ™_a/ and x TOjv-matrices with entries in 

Reduce the rows of the matrix j HiVjSi ■ tj — iAmM modulo a Grobner basis 
for D ■ {Pi, . . . , Pa} C D™" . Force this reduction to be zero by setting the 
coefficients (which are inhomogeneous bilinear polynomials in /i^, i^j) of every 
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standard monomial in every entry to be zero. Collect these relations in the 
ideal Im C K[^ii, . . . , /i^j^j , z/i , . . . , 

4. Similarly, reduce the rows of the matrix "Ylii j ' ■Si ~ idmjv modulo a 
Grobner basis for D ■ {Qi, . . . , Qb\ C Z?™". Force this reduction to be zero 
by setting the coefficients of every standard monomial in every entry to be 
zero, and collect these relations in the ideal In C K[^,i>\. 

5. Put /(y, W) ^ Im + In C K[^l, u]. If I{V, W) contains a unit, return "No" 
and exit. 

6. Otherwise compute an isomorphism X^I^i ^i^i Hom£)(M, N) by finding the 
first T coordinates of any point in the zero locus of I{V, W). For instance, we 
can do this by inductively finding ki ^ K for each i from 1 to r such that 
/(y, W) + {fii — fci, . . . , /ii — fci) is a proper ideal. At each step i, this can 
be accomplished by trying different numbers for ki until a suitable choice is 
found. 

7. Return "Yes" and the isomorphism {J2i=i ^iSi) ■ M N. 



Remark 7.3. Algorithm 7.2 can also be modified to detect whether M is a direct 
summand of N. Namely M is a direct summand of N if and only if the ideal Im 
of step 3 is not the unit ideal. Similarly is a direct summand of M if and only 
if the ideal In of step 4 is not the unit ideal. 

Proof. Reduction of the generic matrix ■ HifjSi-tj—idmM modulo D-{Pi, . . . , Pa} 
in step 3 leads to a generic remainder which depends on the parameters fii , Vj . 
Moreover, since a Grobner basis of D ■ {Pi, . . . ,Pa} is parameter- free, this generic 
remainder has the property that its specialization to a fixed choice of parameters 
fii = Ui, Uj = bj gives the remainder of J2i j o,ibjSi-tj—idmM modulo D-{Pi, . . . , Pa}. 
Thus setting the remainder to zero in step 3 corresponds to deriving conditions 
on the parameters which makes the endomorphism given by j Mi^'jSi • tj 

equal to the identity on M. This is possible if and only if M is a direct sum- 
mand of A^. The analogous statement holds for reduction of j Mi^'j^j • Sj — idmjv 
modulo D ■ {Qi, . . . , Qb} and setting its resulting remainder to zero. Here, setting 
a remainder to zero is equivalent to the vanishing of the coefficients of its stan- 
dard monomials, and we collect these vanishing conditions in the ideal I{V, W) of 

Now a linear combination a.iSi : M N is an isomorphism with inverse 
X bjtj : N M if and only if the composition j a.ibjSi ■ tj is congruent to id^M 
modulo D ■ {Pi, . . . , Pa} and the opposite composition j aibjtj ■ Si is congruent 
to id„i„ modulo D ■ {Qi, . . . , Qb}- Thus the common zeroes (oi, . . . , Ur, bi, . . . , br) 
of I{V,W) correspond to isomorphisms '^^UiSi and their inverses '^jbjtj. In 
particular, if I(V, W) is the entire ring, which we detect by searching for 1 in a 
Grobner basis of I{V, W), then there are no isomorphisms. 

On the other hand if I{V, W) is proper, then M and A^ are isomorphic and we 
obtain an explicit isomorphism from finding any common solution of /(V, W). By 



Lemma 7.1, the invertible homomorphisms from M to A^ are Zariski dense in the 
vector space Homij(M, A^). Hence, a common solution can be explicitly found by by 
intersecting the zero locus of I{V, W) with a suitable number of generic hyperplanes 
{fi = ki}. Because of denseness, each of these hyperplanes can be found in a finite 
number of steps. In other words, if I{y, W) + (^i — /ci , . . . , /i — fci_i) is proper, then 
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there are only finitely many ki for which the sum I{V, W) + {fii — fci, . . . , /x — fcj) is 
the unit ideal. □ 

Remark 7.4. Once we have specialized the /Xj in a common solution of I{V, W), 
then the Vj are determined because of the bilinear nature of the relations (which 
gives linear relations for the i/j once all fXi are chosen). This also means that if 
there is any solution, then the /Uj are rational functions in the vj and vice versa. In 
particular, if (j> £ Hom£)(Af, N) is defined over the field K then (j>~^ is defined over 
K as well and no field extensions are required. We now give two simple examples, 
one where M and N are isomorphic, and one where they are not. 

Example 7.5. Let n = 1 and M = iV = D/D ■ d^. One checks that V = W = 
Homij (M,iV) is generated by the 4 morphisms si = -(9), S2 = -{xd), S3 = •(!), 
and S4 = ■{x'^d — x). We obtain the generic morphism 

4 4 

i^i'^jtj -Si-i = (/i3i/3 - /xii^4 - 1) 

+ (-^4^3 - M2J^4 - M3Z^4)a; 
+ (/"3J^l +Ml!^2 + ^1^^3)9 

+ (-^4i^l + M2J^2 + M3i^2 + + jJLlVijxd 

plus 9 other terms which are va. D ■ independently of the parameters. 

Hence in order for Y^l^^ iJ^iSi to be an isomorphism, the yUj need to be part of a 
solution to the ideal 

I{V, W) = (/Usfs - fj.1i/4 - 1, 

-^4J^3 - M2i^4 - M3«^4, 
/i3J^l + AilJ^2 + Ml^^3, 

-1^41^1 + fl2l^2 + H3^2 + 1^21^3 + l^l^i, 
/U4f3 +M2J^4 + At3!^4)- 

This ideal is not the unit ideal and has degree 8. Hence there are isomorphisms 

between M and iV. Pick "at random" jii = 1, /i2 = 2, and /Z3 = 0. Then the ideal 
I{V, VF) + (jLti - 1, /X2 - 2, H3 - 0) equals the ideal (/^i - 1,^2 - 2, /X3, 1^4 + 1, 1^2 + 
2^35 ^^1 + 5Z^3,At4^'3 — 2). We see that we have to avoid fj.4 = but otherwise have 
complete choice. 

Example 7.6. Let n = 1, M = D/D ■ d"^, and N = D/D ■ d. One checks that V = 

Homi3(iV, M) is generated by h = -{d) and -{xd-l) while W = Homz3(Af, N) 
is generated by si = •(!) and S2 = -(x). The sum ^ fJ^i^jSi ■ tj takes the form 

fl2'/2X^d + {filV2 + ^l2Vl)xd + flll^ld - {pil'2 + A*2^'2). 

Modulo D ■ d we want this to be 1, so wc get the relation 

M2J^i - ^J'll^2 = 1- 

We note that this equation has plenty of solutions, which means that M can be 
realized as a summand of A''. On the other hand, the sum ^ HiUjij ■ Si takes the 
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form 

Modulo D ■ we want this to be 1, so we get the relations 

-^li^2 = 1, 

fiiui = 0, 

^lJ/2+M2t^l = 0, 
/i2i^2 = 0. 

Putting all the equations together, we obtain the unit ideal, and hence M and N 
are not isomorphic. 

For M and N isomorphic, we now give a method to find all possible isomor- 
phisms, that is we study the units Isod{M) in the endomorphism ring V = End£)(M). 

Lemma 7.7. The isomorphism set lso£i{M) is a smooth affine variety which is 
connected if K is algebraically closed. 

Proof. As we have seen, 1sod{M) is isomorphic to the nonempty affine variety 
Var(M) = Var(/(V, V^)) defined in the variables Hi^Vj. Here, a point of Var(M) 
with coordinates (/ii, . . . , /i^, z^i, . . . , Vt) corresponds to the isomorphism ^11=1 f^i^i- 
Now any isomorphism : M — > M induces an isomorphism from the variety to 
itself, sending {ij.,v) to {n' ,v') where ^HiSi o (jj = J2l^i^i- This action is regular 
in (j) (since we showed /i is rational in i/), and transitive since ip G Hom^i (M,M) 
equals (-0 ° o It follows that Var(M) is a smooth variety because it is a 
homogeneous space over itself via a transitive action. As we have seen in Lemma 



7.1, Iso(M) is Zariski open in EndD(M) (which is an affine space and therefore 



normal) and hence connected if K is algebraically closed. □ 

Since the isomorphisms Iso d{M) arc Zariski open in EndD(M), one can ask 
for a method to compute the equations defining the complementary closed set of 
non-isomorphisms . 

Definition 7.8. The ideal in Klv] that determines the closed set End£)(Af) \ 
Isod(M) of non-isomorphisms of M is called A(A/), the defect ideal. 

Algorithm 7.9. (Computing the defect ideal) 
Input: Generators for a holonomic D-modvle M. 

Output: The defect ideal A(M) defining the non-isomorphisms of End£)(M). 

1. Perform Steps 1 through 4 of Algorithm with M — N as input to obtain 
the ideal I{V, V) C K[{fi,} U {i^^}]. 

2. Regard each of the C generators of I{V, V) as a linear inhomogeneous equation 
in the variables fii with coefficients involving Vj as parameters, and collect all 
these equations in a single matrix equation A ■ ji = b, A G K[i']''^'^ . 

3. Compute all r x t minors of A and collect them in an ideal A(M) C Kli']. 

4. Return A(M). 

Proof. A point v corresponds to an isomorphism with inverse fi if and only if the 
system A-jj. ^ b has exactly one solution for /x. This is equivalent to the C x r matrix 
A having rank t and the augmented matrix {A\b) also having rank r. The matrix 
A will have rank r if and only if any one of its r x r minors is nonzero. Similarly, 
the augmented matrix {A\b) will also have rank t if in addition all (r + 1) x (r + 1) 
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minors vanish. We claim that each (r+ 1) x (r + 1) minor of {A\b) must actually be 
identically zero. Otherwise it would impose an algebraic condition which must be 
satisfied by the isomorphisms in the coordinates v of Endc(Af). But this cannot 
happen since the isomorphisms are an open set by Lemma 7.1. Thus we have shown 
that the space of non-isomorphisms v is defined by the equations obtained from the 
vanishing of all r x t minors of ^. □ 



Remark 7.10. In Lemma 7.1, we saw that modulo the Jacobson radical, which is 
a linear subspace, then EndD(M) is the product of simple iC- algebras. Moreover 
when K is algebraically closed, then a simple X-algebra is a matrix algebra. It 
follows that in an algebraic closure of K, the radical of A(M) is generated by linear 
forms corresponding to the Jacobson radical and a single determinant which is the 
product of the determinants of the matrix algebras. However, we have not yet 
understood what happens when the field of input K is not algebraically closed. 



Optimistically, we hope that Algorithm 7.£ produces an ideal A(M) whose radical 
over K is also generated by linear forms and a single determinant. The following 
example shows at least that A{AI) might not be radical. 



Example 7.11. Let us look at our Example 7X. There the system of equations 

{-^11^4 + ^^3*^3 = l,-At2J/4-A«3l'4-M4l^3 = 0, fJ'liv2 + V3) + ^i^l^i = 0, /iiZ/4 + A*2 (^^2 + 

^3) + fJ-3i^2 ~ fJ-ii^i = 0} can be rewritten in the form A ■ fj, — b as, 
I 



-V4 



V2 + 1^3 



^3 
1^2 



'V3 
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( ,1] 








tJ-2 









M3 
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\ / 




^0/ 



In order to assure that A has full rank we need the determinant 



A{M) = v^v^ + 2v2vl 



2V1V2V3VA 



21^1 1/3 2^4 + ViV^ 



to be nonzero. We conclude that the locus of not invertible morphisms is given by 
the vanishing of the determinant of A(M). Note also that A(M) — {v2V^ + '^l + 
1/11/4)^ which in particular is not square free. 

We end by discussing the endomorphism ring E = Endi5(Ai^). 

Remark 7.12. A well-known application of i? = EndD(Af) is towards decompo- 
sitions of M . By the KruU-Schmidt-Azumaya theorem, a Z?- module M has a de- 
composition into a direct sum of indecomposable submodules (meaning that they 
cannot be further decomposed into a direct sum of nonzero submodules) , and any 
such decomposition is unique up to re-ordering and isomorphism (see e.g. [|[ The- 
orem 19.21]). There is a bijective correspondence between (1) the decompositions 
of M into a direct sum of submodules and (2) the decompositions of the identity 
element 1 ^ ei + ■ ■ ■ + eg oi E into pairwise orthogonal idempotents ^, Theorem 
1.7.2]. The correspondence is gotten by taking a set of orthogonal idempotents 
{ei, . . . , e^} and producing the decomposition M = ei • M © • • • ® • M. Thus, an 
algorithm which produces a full set of orthogonal idempotents for the X-algebra 
End£)(Af) combined with Algorithm 5.1 would give a method to decompose holo- 
nomic D-modules into indecomposables. 

Computation in finite-dimensional iiT-algebras E has recently been an area of 
active research. When K is a, number field, early work of Friedl and Ronyai provides 
polynomial-time algorithms to decompose E into simple algebras if E is semi-simple 
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and to find the radical of E in general Q. When K — C, Eberly has given Las 
Vegas polynomial time algorithms to find the decomposition of a simple algebra as 
a full matrix ring |4| . We should also mention that the radical of E is independent 
of field extension of K while the decomposition into simple algebras depends upon 
the field K. Thus if we are willing to use K — C, then a full set of orthogonal 
idempotents can indeed be algorithmically computed. 

Let us also describe another method based on computational algebraic geometry 
to obtain information on the invariants di in the decomposition 

d 

E/ Jsic{E) (^kK = Y[ EndjiiK'^^) 
1 

where K denotes the algebraic closure of K. We will compute the de Rham coho- 
mology groups of the complement of Var(A(M)) in EndD(M) = C"^ = Spec(C[j/]) 
using the algorithm developed by the second author in |l6| . This algorithm will in 
some sense allow us to pretend that K is already algebraically closed. Namely, the 
algorithm can be used on input defined over any computable subfield of the com- 
plex numbers, but always computes dimc{H*j^{C" \ Y, C)). What we now need is a 
method that enables us to sort out the di from the Betti numbers of C'^\Var(A(M)). 

Consider E ®k C. As we have shown, its units are (homotopy equivalent to) 
the units of a product of the form JliLi Endc(C'*'). The non- isomorphisms in 
each factor are given by the vanishing of the appropriate determinant, and the 
isomorphisms are just the elements of the general linear group Gl{di, C). 

The cohomology of Gl{n, C) is well understood and best expressed for our pur- 
poses in terms of the Poincare polynomial. 

Definition 7.13. Let ho, . . . ,hi,... be the dimensions of the de Rham cohomology 
groups of a complex manifold T. Then the Poincare series (polynomial) Pxiq) is 
defined by 

i>Q 

The Poincare polynomial behaves very nicely under products M2 = Mi x A/3 of 
manifolds: 

Pmi (9) • -Pa/3 (9) = Pm2 (?)• 
An old result (||l9|, Theorems T.ff.A and 8.16.B) states that PGi(rf..C)(9) = 
n^=i(l + t^-'^^). Hence the Poincare polynomial of a product of general linear 
groups lY^^^ PGi(d,,C){(l) equals 

d di 

(If) n 11(1 + ^ _^ t)'^^^>o 1 . (1 + t3)E<i.>i 1 . (1 _^ t5)E.,>2 1 . , _ 

Thus in order to compute the di one then has the following algorithm. 
Algorithm 7.14. 

Input: Generators and relations for the left module M. 
Output: The invariants di associated to End£)(M). 

1. Compute the defect ideal A(M) C Klv] by using Algorithm |7.9| . 

2. Compute the dimensions hk = dimc(i?dfl(C" \ Var(A(Af))). 
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3. Factor the polynomial Pison{i^i)il) ■— J2 hk(l^ into 

+ + .....(l + g2i-l)fe,_ 

4. Compute the di by comparing the expansion from the previous item with 
equation (O). Return the di. 



Example 7.15. Continuing our Example 7.11 we compute the de Rham cohomol 



ogy groups of the complement oiYa,T:{v2V^ + 2v2V^ + v'^ + 2viV2V3,V4: + 2vivi^Vi + vivj) 
in C*. Using Macaulay2 one obtains ~ hi = = = 1 and all other hk van- 
ish. Then the Poincare polynomial is 1 + q + + q'^ = (1 + + q^). This means 
that d — 1 and di = 2. 

8. Appendix 

In this section, we provide a short survey of the ideas that lead to an algorithm 
for restriction and then, mostly for purposes of reference here and otherwise, list 
an algorithm to compute integration. All the main ideas are taken from 0, |ri|, ^ 

Definitions 8.1. Fix an integer d with < d < n and set H = Var(a;i, . . . , Xd)- 
For a G Z", we set uh = (ai, . . . , ad, 0, . . . ,0). 

On the ring D we define the Vd-filtration F^{D) as the X-linear span of all 
operators a;" 9^ for which |a/f| + fc > \(3h\- More generally, on a free D-module 
A = ®*=i-D ■ Cj we define 

where m is an element of Z™ . We shall call m the shift vector. A shift vector is tied 
to a fixed set of generators. 

We define the Vd-degree of an operator P G ^[m], Vd deg(P[m]), to be the smallest 
k such that P e F^(^[m]). 

If M is a quotient of the free D-module A — ®\D ■ ej, M = A/ 1, we define 
the Vd-filtration on M by F^(Af [m]) = F^{A[m]) + I. For submodules iV of A we 
define the Vd-filtration by intersection: F^{N[m]) = F^{A[m]) D N. 

If M is a submodule of the free module ^[m], then a Vd-strict Grobner basis 
or a Vd-Grobner basis for AI is a set of generators {mi,... ,mr} for M which 
satisfies: for all to G Af we can find G D such that to = ^aiTO^ and 

Vd deg(Q;iTOi[tn]) < Vd deg(TO[m]) for all i. 

Definitions 8.2. A complex of free Z3-modules • • • — > A^~^ '^-^ A'- ^ A^'^^ — > • • • 
is said to be Vd-adapted at A^ with respect to certain shift vectors mi_i, mi, rrii+i 
if 

0^ (F^(A^[m,]))cF^(^'+i[m.+i]) 

and also 
for all k. 

We shall say that the complex is Vd-strict at if it is Vd-adapted at and 
moreover 

im((/)^-i) ni^^(yl'[m.]) = im(0»-i|^.(^.-,[^^_^])) 
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for all k. 

For 1 < d < n we set 6d = xidi + . . . + Xddd and = 0. Recall that a D-module 
M[m] = A[m]/I is called specializable to H if there is a polynomial b{s) in a single 
variable such that 

(12) b{ed + k)-F^{M[m]) C F^-\M[xa\) 

for all k (cf. ^^)- Introducing 

gr^(M[m]) = (F^(M[m]))/(F^-i(Af[m])), 
this can be written as 

b{ed + k)-gy\j{M[m])^Q. 



The monic polynomial b{9) of least degree satisfying an equation of the type (12) 
is called the b-function for restriction of M[m] to H. 

Remark 8.3. Spccializability descends to quotients and submodules. Namely, 
assume that A/[m] = (A//)[m] is specializable and A^[m] = (A'//)[tn] is a sub- 
module of M (where I A' C A). Let b{s) be a polynomial that satisfies 
b{0d + k) ■ F^{A[m]) C F^-^{A[m]) + /. Then clearly 6(0^ + k) ■ F^{A[m\) C 
F^~^(A[m]) + A' as well and hence {M/N)[m] is specializable to H. On the other 
hand, if P' £ F^{A'[m]) = F^{A[m]) n A' then b{e + k) ■ P ^ Q + Q' where 
Q e i^^"i(A[m]) and Q' e I and hence Q e F^-\A[m\) n A' ^ F^{A'[m]). This 
implies that N is also specializable and we see that the 6-functions for restriction 
of A^[m] and for {M/N)[m] divide the 6-function for restriction of M[m\ to H . 

Notice that independently of d, gr^(Z?[0]) = _D, as a ring. 

It has been shown by Oaku and Takayama in |l^ (Proposition 3.8 and following 
remarks) how to compute Vd-strict Grobner bases, and for any D-module M posi- 
tioned in degree b a free V^-strict resolution {A' [m,], 4>') of Af [m;,], — (B\^D, ri = 
if i > 6. The construction given in [l0[| allows for arbitrary m;,. 

The method employed is to construct a free resolution with the usual technique 
of finding a Grobner basis for ker(yl* A*"*"^) and calculating the syzygies on this 
basis. The trick is to impose an order that refines the partial ordering given by 
Vd-degree, together with a homogcnization technique. 

In jl^ was given an algorithm to compute y„-strict resolutions for right bounded 
complexes, based on Eilenberg-MacLane resolutions. In [|l^ an improved algorithm 
is given that usually computes a much smaller resolution and is also easier to 
implement. We shall assume that the reader is familiar with the techniques from 

(io[ [111110. 



An idea first stated in and further developed in |T^] yields a theorem which 
in order to state we need to introduce some more terminology for. 

Definition 8.4. Let = D/{xi, ... , a;^} ■ D and fi^ = D/{di, ... ,dd}- D. 

The restriction of the complex A*[m,] to the subspace H is the complex tld 
A*[Tn,] considered as a complex in the category of i^(a;c(+i, . . . , x„, 9„)-modules 

The integration of A*[m,] along H is the complex ild ®|3 A*[tTi,] considered as 
a complex in that same category. 

We need to make a convention about the V^-filtration on tensor products over 

D. 
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Definition 8.5. If ^[m] is a free i7-graded 13- module with shift vector m then 
Cld (X)D A[m\ is filtered by F^{Cld 0_d A[m]) := the /f-span of {P(g)zj Q\Vddeg{P) + 
Vddeg{Q[m]) < k}. Note that as f2<j equals K[di, .. . ,dd]{xd+i,dd+i, ■ ■ ■ ,x„,a„)as 
right _D-module, F^{rid'S)D A[m]) equals the free K{xd+i,dd+i, ■ ■ • , Xn, 9„)-module 
on the symbols {(Pi, . . . , PrkniA))\Pj e K[di, ■■■ , dd],degg{Pj) <k- m{j) Vj}. 

If A*[m,] is a Vd-strict complex, we denote by F^{(l,d 65 ^•[m,]) the complex 
whose modules are the P^(f2d(X)A*[mi]) as defined above, and the maps are induced 
from A*. 



Algorithm 8.6. Let (A* [m,], (5*) be a V^-strict complex of free D-modules. The 
restriction of A* [m,] to H = Var(a;i, . . . , Xd), interpreted as a complex of modules 
over K{xd+i, dd+i, ■ ■ ■ , a;„, 9„), can be computed as follows: 

1. Compute the fe- function fo^'lm.] (*) for restriction of A*[m,] to H. 

2. Find integers fco, fci with (&A'[m.] (fc) — 0, fc G Z) ^ (fco < k < ki). 

3. (id f^j^ A' is quasi- isomorphic to the complex 



(13) 



This is a complex of free finitely generated K{xd+i, dd+i, • • • , a;„, 9„)-modules. □ 

Let us now consider the question how to compute fid <E)j^ C . This problem is 
intimately related to the restriction algorithm. The reason is the Fourier transform^ 
which is an algebra automorphism from D to itself and defined as follows: 

Td{xi) = di, T{di) = -x^. 

The perhaps surprising minus sign is required to keep the Leibniz relation 0:^9^-1-1 — 
diXi intact. The Fourier transform can be used to define an equivalence of the 
category of left D-modules with itself via "extension of scalars" : Td{M) := D®dM 
where D is on the left considered as just D while on the right D acts on D via Td- 
So for example if 1 (8) to G J-d{M) then Xi-l®m = Xi®m=l® {—di ■ m). 

If we apply J^d to the integration problem we are reduced to computing the 
restriction of TdiC') since J-d{^d) = (id- We are led to introduce therefore a 
Vd-filtration which is defined by 

F^H{D)^Td{F^H{D))- 
This extends just as the Vd-filtration to submodules and quotients of shifted free 
modules. One also defines a b-function for integration of the complex C", be, as 
the 6- function for restriction of the complex Td{C*). Let us illustrate this concept 
with an 

Example 8.7. Suppose n — 2, M = D2/D2 ■ {81,82) = i^[a;i,a::2] and d = 1. 
The 6- function 6(6*1) for restriction of the complex X* with X'^ = M and X* = 
otherwise is &(6'i) = 6*1 because xidi ■ F%{D2) C F^\D2) + D2 ■ {91,52}. On the 
other hand, the 6-function b{Oi) for integration is 6*1 + 1 because (xiC?i + l)-F^(-D2) ^ 

F^\D2)+D2-{xi,X2}. 



Theorem 3.6 implies then the following algorithm. 



Algorithm 8.8. Let (A'[m,],(5*) be a V^-strict complex of free D-modules. The 
integration of A'[m,] along di, . . . ,dd, interpreted as a complex of modules over 
K{xd+i, dd+i, ■ ■ ■ , Xn, dn), can be computed as follows: 
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1. Compute the 6- function ^A'lm.lls) for integration of A*[m,] along H. 

2. Find integers fco, fci with (&A'[m.](fc) = 0, fc e Z) => (fco < fc < fci). 

3. Vid ®^ A* is quasi-isomorphic to the complex 

This is a complex of free finitely generated K{xd+i, dd+i, ■ ■ ■ ,Xn, 9„)-modules. 
Example 8.9. K[xi,X2] has a Vi-strict free resolution 

D[-i] i^ei?[o,-i] IzM D[Q]. 

Continuing our example, we want to find fci, satisfying the condition ((/s + 1) = 
0, fc G Z) ^ /co < < ^1- Clearly — ki — —\ should be chosen. 

Then the integration of M along di is, according to the theorem, quasi-isomorphic 
to the complex 







Fh\D2) Fh\D2) Fh\D2) 

Since Fj^^{D2) is the span of all monomials of D2 with positive Vi-degree and 
F^^{D2) is spanned by those of Vi-degree exceeding 1, the complex above is (with 
Di=K{x2,d2)) 

>L>i-l^L>i-l©0^0^0^--- 

the cohomology of which is exactly K[x2], shifted cohomologically by one relative 
to the input. 
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